DIAMETERS OF RANDOM CIRCULANT GRAPHS 



JENS MARKLOF AND ANDREAS STROMBERGSSON 

Abstract. The diameter of a graph measures the maximal distance between any pair of 
vertices. The diameters of many small-world networks, as well as a variety of other random 
graph models, grow logarithmically in the number of nodes. In contrast, the worst connected 
networks are cycles whose diameters increase linearly in the number of nodes. In the present 
study we consider an intermediate class of examples: Cayley graphs of cyclic groups, also 
known as circulant graphs or multi-loop networks. We show that the diameter of a random 
circulant 2fc-regular graph with n vertices scales as n^^'^ , and establish a limit theorem for 
the distribution of their diameters. We obtain analogous results for the distribution of the 
average distance and higher moments. 



1. Introduction 

The diameter of a graph is the largest distance between any pair of vertices, and is a 
popular measure for the connectedness of a network. Many models of small- world networks, 
for example, have diameters that grow slowly (i.e., logarithmically) with the total number of 
nodes [9], [19]. The same phenomenon is observed for a wide variety of other random graph 
models, and has been proved rigorously in many instances [6], [71, [12], [15], [18], [27], [33] . 
The worst connected networks are cycles, whose diameters increase linearly with the number 
of vertices. Here, connectedness is dramatically improved by additionally linking every vertex 
with a random partner; the logarithmic growth of the diameter is then recovered [8]. 

In the present paper we consider a more regular generalization, the circulant graphs (often 
also called multi-loop networks) which comprise an interwoven assembly of cycles (Figs. [l| [2] 
left). We will show that the diameter of a random 2/c-regular circulant graph with n vertices 
scales as n^^^ , and prove a limit theorem for the distribution of diameters of such graphs; the 
existence of a limit distribution was recently conjectured in |2j. Analogous results hold for 
the distribution of the average distance in a circulant graph and related quantities, see Sec. 
[5] for details. It is interesting to note that an algebraic scaling of the diameter has also been 
observed for the largest connected component of the critical Erdos-Renyi random graph |31j ; 
here the scaling factor is n^/^. 

We furthermore establish corresponding results for circulant digraphs (cf. Figs. [T| [2] right), 
where the limit distribution of diameters turns out to coincide with the limit distribution 
of Frobenius numbers in d = k + 1 variables studied in |28j. The connection of these two 
objects has been exploited previously [3], [32], [35], [39]. As for the Frobenius problem [23] . 
the question of calculating the diameter of circulant graphs can be transformed to a problem 
in the geometry of numbers [11], [41]. We will use a particularly transparent approach that 
identifies circulant graphs with lattice graphs on flat tori [13j, [16], and then employ the 
ergodic-theoretic method developed in [28] to prove the existence of the limit distribution of 
diameters. 
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Figure 1. The 4-regular circulant graph C8(2,3) and the chxulant digraphs 
C^(2,3), Cg{2,5). The corresponding diameters are 2, 3 and 4, respectively. 




Figure 2. The 3-regular circulant graph Cio(2,5) and the circulant digraphs 
C^(2,5), CfQ{5,8). The corresponding diameters are 3, 5 and 5, respectively. 

Let us fix an integer vector a = (ai, . . . ,ak) with distinct positive coefficients < ai < . . . < 
Qfc < f • We construct a graph Cn{a) with n vertices 0, 1, 2, . . . , n — 1, by connecting vertex i 
and j whenever |i — j| = ah mod n for some h S {1, . . . , fc}. Because the adjacency matrix of 
this graph is circulant, Cn{a) is called a circulant graph. If Ofc < §, then C„(a) is 2A;-regular, 
i.e., every vertex has precisely 2k neighbours. If = then Cn{a) is {2k — l)-regular. It is 
easy to see that Cn{o) is connected if and only if gcd(ai, . . . , a^, n) = 1. In this case Cn{o,) is 
the (undirected) Cayley graph of 'L/n'L with respect to the generating set {±ai, . . . , ^a^}- 

To construct a directed circulant graph [circulant digraph for short) choose an integer vector 
a = (ai,...,afc) with distinct positive coefficients < ai < . . . < < n. The circulant 
digraph C^{a) is defined to have an edge from i to j whenever j — i = mod n for some 
/i G {1, . . . , k}. In C^{a), every vertex has precisely k outgoing and k incoming edges. C^{a) 
is strongly connected if and only if gcd(ai, . . . ,ak,n) = 1. In this case C+(a) is the directed 
Cayley graph of Z/nZ with respect to the generating set {ai, . . . , Ofc}. 

Fix a vector £ = [£i, . . . , i^) S M^q. We endow our circulant (di-)graph with a (quasi-)metric 
by stipulating that the edge from i to j = i + mod n has length i^. We denote the 
corresponding metric graphs by C„(£, a) and C^(£, a), respectively. The distance d{i,j) 
between two vertices is the length of the shortest path from i to j. The diameter is the 
maximal distance between any pair of vertices, 

(1.1) diam = maxd(i, j). 

To define an ensemble of random circulant graphs, we set 

(1.2) d+ ■.= {xeR''+^ :0<xi<...<Xk<Xk+i}; ^ := n {xk < ^Xk+i}, 

and then in the directed case we fix an arbitrary bounded subset T> C 5"*" with nonempty 
interior and boundary of Lebesgue measure zero; in the undirected case we fix an arbitrary 
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Figure 3. Distribution of diameters n^^/'^(diam C^(€, a) + e- £) for circulant 
digraphs with k = 2 and i = e := (1, 1) vs. Ustinov's distribution P2{R) in 
(1.12). The numerical computations assume T) = f^{x^ < 1} and T = 1000. 



bounded subset V <Z subject to the same conditions. Denote by N the set of integer 
with positive coprime coefficients (i.e., the greatest common divisor of all 



vectors in 



sfc+l 



coefficients is one). The numbers (ot,n) E fi^^^ defining Cn{£,o,) or C^{£,a) are then picked 
uniformly at random from the dilated set TT? (T > 0). Note here that N'^"'"^ nTP is nonempty 
for all large T; in fact 

vol(P) 



(1.3) 



#{(a,n) G n''+^ nrp} 



as T 



oo. 



C{k + l) 

Our first main theorem shows that the (properly scaled) diameter of a random circulant 
digraph has a limit distribution which is independent of the choice of P. In order to describe 
this limit distribution, we introduce some further notation. For a given closed bounded convex 
set K of nonzero volume in M.^ and a (fc-dimensional) lattice L C M^, we denote by p{K,L) 
the covering radius of K with respect to L, i.e. the smallest positive real number r such that 
the translates of rK by the vectors of L cover all of M'^: 

(1.4) p{K,L) = mi{r>0 : rK + L = R''}. 

Let Xj. be the set of all lattices L C M'^ of covolume one, and let /io be the SL(/c,M) invariant 
probability measure on X^. Also let A be the simplex 

(1.5) A = {x£R'^Q:xi + ... + Xk<l}. 

Theorem 1. Let k > 2. Then for any I G M^q and any bounded set D C 5^ with nonempty 
interior and boundary of Lebesgue measure zero, we have convergence in distribution 

diamC+(£, a) 



(1.6) 



/o(A, L) as T ^ oo, 



(n£i---4)i/'= 

where the random variable in the left-hand side is defined by taking (a, n) uniformly at random 
in N'^'^^nTP, and the random variable in the right-hand side is defined by taking L at random 
in Xk according to fiQ. 



Remark 1.1. The limit distribution in Theorem [T] is the same as the limit distribution for 
Frobenius numbers in d = A; + 1 variables found in ^28j, and our proof depends crucially 
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Figure 4. Distribution of diameters n ^^'^ diam C„(e, a) for circulant graphs 
with A; = 2 vs. our formula (1.20). The numerical computations assume T> = 
Sn{x3< 1} and T = 1000. 



on the equidistribution result proved in [28, Thms. 6, 7]. Let Pk{R) be the complementary 
distribution function of p{A,L), viz. 

(1.7) Pk{R) := fio{{L e Xk : p{A,L)>R}). 

{Pk{R) = ^d{R) in the notation of [28j.) It was proved in [25 that Pk{R) is continuous for 
any fixed k > 2. Hence, recalling also (1.3), the statement of Theorem [l] is equivalent with 
the statement that for any i? > we have 



lim --^#Ua,n) nTV 



diamC^(£, a) 



> R 



vol(P) 

0^+1) 



PkiR). 



We also remark that Li [26J has recently proved effective versions of the equidistribution 
results in [5H]. Using Li's work it should be possible to also prove effective versions of our 
Theorems [T| [2j as well as Theorems [3| |4] in Section [2| 

Remark 1.2. In analogy with the case of Frobenius numbers [1], we also obtain the following 
sharp lower bound, writing e := (1,...,1) sM'^, 



(1.9) 



diamC+(.£, a) + e-i 



> pk, with pk := inf p{A,L). 



It follows from the description in Remark |1.1| that 

(1.10) Pk{R) = 1 for 0<R<pk, and < PkiR) < 1 for R > pk. 

It is proved in [1] that pk > (kl)^^^, and in fact for k large, pk is not much larger than (/c!)^/'^; 
indeed pk < (kl)^ ^ {1 + 0{k-^ log k)) (cf. [Ml Sec. 9], [20], [Ml). Also for k large, the limit 
distribution described by Pk{R) has almost all of its mass concentrated between {k\)^^^ and 
1.757 • (/c!)i/^. In fact, for any fixed a > 1 + rjo, where rjo = 0.756 ... is the unique real root of 
elogr/ + T] = 0, Pfc(a(A;!)i/'=) tends to zero with an exponential rate as — )• oo \38\ Thm. 4.1]. 
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Remark 1.3. For k fixed and R large, 

(1-11) Pk{R) = ^^R^' + Ok{R-'-'-^^). 

This asymptotic formula is proved in [38| Thm. 1.2]. The upper bound Pk{R) < -R"^ had 
previously been proved in 



Remark 1.4. For k = 2, Theorem [T] has been proved by Ustinov by different methods, see 
the last section of [39]. This paper also computes an explicit formula for the limit density 
Pk{R) = ~'^Pk{R) (which coincides with the distribution of Frobenius numbers for three 
variables): 

'0 (0<i?<\/3) 

(1.12) p2{R) = < 



i^(^_V4^) iV3<R<2) 
i|(i?V3arccos(^?i5^) + iVir^log(f5|)) (i? > 2). 



We give an alternative proof of this formula, deriving it as a consequence of ( |1.7[ ) , in Section 
lU below. 

We now turn to the case of undirected circulant graphs. The following theorem says in 
particular that, as in the directed case, the limit distribution for the diameter is independent 
of the choice of V. Let *P be the (regular) polytope 

(1.13) ^ = {x eR^ : \xi\+ . . . + \xk\ <l}. 

This is a /c-dimensional cross-polytope, cf. |14| : in particular ^ is a square for k = 2 and an 
octahedron for k = 3. 

Theorem 2. Let k > 2. Then for any i G M^q and any bounded set P C 5 with nonempty 
interior and boundary of Lebesgue measure zero, we have convergence in distribution 

(1.14) (n£,---4)V^ ^^^^''-^ asT^oo, 

where the random variable in the left-hand side is defined by taking (a, n) uniformly at random 
in N^^^ nTV , and the random variable in the right-hand side is defined by taking L at random 
in Xk according to fiQ. 

Remark 1.5. Let Pfc(-R) be the complementary distribution function of p{^,L), viz. 

(1.15) Pk{R) := ^o{{L e Xk : p{%L)>R}). 



This function is continuous (cf. Section [3.1| below), and hence, recalling also (1.3), the state^ 



ment of Theorem [2] is equivalent with the statement that for any i? > we have 

(1.10, ^*{K„,.S-n™ : f=^>4^^n,.). 

Remark 1.6. We have the lower bound (cf. Proposition [l] and Lemma [4] below) 

diam Cr,(£, a) + ie • £ 
(1-17) , , ,ul > Pk, with pk := inf p(q3, L). 

(Recall e := (1,...,1) G M'^). Also the distribution described by Pk{R) has support exactly 
in the interval [pk,oo), in analogy with (1.10). Since any covering of M.^ has density at least 



one we have 

(1.18) Pfc>vol(q3)-Vfc= i(fc!)iA. 

In fact (1.18) holds with equality for k = 2; p2 = since there exist lattice coverings of ^ 



»2 



by squares without any overlap; however for every A; > 3 we have strict inequality in (1.18); 
cf. Section [3^ below. We also have pk < l(/c!)i/''(l + 0{k-^logk)) (again cf. [M Sec. 9], 
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Figure 5. Distribution of the shortest cycle length n~^^'' sclC^ (e, a) for cir- 
culant digraphs with k = 2 vs. the probability density P2,sci(-R) discussed in 
Section [s] The numerical computations assume P = n {xs < 1} and 
T = 1000. 

|20| . |34j). and for k large, the limit distribution described by Pk{R) has almost all of its 
mass concentrated between ^(kl)^^'' and 1.757 • ^(A;!)^/*'', in the same sense as for Pk{R) [38^ 
Thm. 4.1]. 



Remark 1.7. For k fixed and R large, we will show in Section 3.3 that 
(1-19) Pk{R) = ^)+Ok{R-'-'-^^). 

Remark 1.8. In the case k = 2, the limit density pk{R) = ~ZRPk{R) can be calculated 
explicity; we will show in Section [4] that 



(1.20) P2{R) 




(0 < < ^) 

^#iog(^)) iR>^)- 



The outline of the paper is as follows. In Section [2] we prove Theorems [T] and [2j by realizing 
the circulant graphs as lattice graphs on flat tori, and applying the central equidistribution 
result proved in |28j . In Section [3] we prove the assertions which we have made about the 
limit distribution in Theorem [2| viz. that the distribution function R i— t- Pk{R) is continuous, 
that we have strict inequality pk > ^(A;!)^/^ for every A; > 3, and that Pk{R) has the precise 
polynomial decay as given by pTTol ). In Section [4] we prove the explicit formula for p2{R), 



and also give a new proof of the explicit formula for p2 (R) ■ Finally in Section [5] we discuss a 
number of natural extensions and variations of Theorems [T] and O 

Acknowledgements. We are grateful to Svante Janson for inspiring and helpful discussions. 

2. Lattice graphs on flat tori and their continuum limit 

In this section we will prove Theorems [T] and [2j The first step is to realize an arbitrary 
circulant graph as a lattice graph on a flat torus. This has previously been used in [16] and 
[13\; we here give an alternative presentation, adapted so as to make the equidistribution 
results from |28j apply in a transparent fashion. 
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Figure 6. Distribution of the shortest non-trivial cycle length 
n~^^'' sclCn{e,a) for circulant graphs with A; = 2 vs. the probability 
density p2,sci(-R) discussed in Section [sj The numerical computations assume 
p = g- n {xs < 1} and r = 1000. 



2.1. Directed lattice graphs. Let LG^ = (Z^,i?) be the standard directed lattice graph 
with vertex set Z'^; the edge set E comprises all directed edges of the form (m, m + eh) where 
m £ Z'^ and ei,...,efc is the standard basis. We define a quasimetric on LG^ by fixing 
£ = {£i, . . . ,£k) £ I^>o assigning length to every edge of the form {m,m + e/^). The 
distance from vertex m to n in LC^ is then given by 



(2.1) d{m,n) 



m) ■ £ if n — m G Z>q, 
otherwise. 



If A is a sublattice of we define the quotient lattice graph LG^ / A as the digraph with 
vertex set Z'^/A and edge set 

(2.2) {(m + A, m + + A) : m G Z'', h = l,...,k}. 

(Note that edges of the form (m + A, m + A) correspond to loops.) The distance from vertex 
m + A to n + A in LG^ /A is 

(2.3) d(,„ + A.n + A) = |'™«"-'" + '^)"^>»)^'^ it(n-m + A)nZ|„^«. 

I oo otherwise. 

Set d = k + 1. Given (a, n) = (oi, . . . , ak, n) G N*^ with < ai < • • • < < n, we introduce 
the following sublattices of U^: 

(2.4) A„ = Z'' X nZ and A„(a) = A„n(a), 
where 

(2.5) u{a):=(^^ ^J") G SL(d, Z). 

For a subset y C M'' we denote by Yq the set Y n (M^ x {0}); we view Yq as a subset of M'^. 

Lemma 1. The set An(a)o is a sublattice of of index n; furthermore the quasimetric 
digraphs LGt / ^n{0')o oind G^{£,a) are isomorphic. 
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Proof. An integer vector m G Z*^ lies in A„(a)o if and only if (m, 0) E A„(a), and this holds 
if and only if m • a = mod n. In other words A„(a)o is the kernel of the homomorphism 
m ^ m ■ a mod n from Z'^ onto TLjuL. Hence A„(a)o is indeed a sublattice of Z'^ of index 
n, and the map just considered induces an isomorphism J : Z'^/A„(a)o — s- Z/nZ. Note that 
J(e/i + An(a)o) = a/j mod n; hence the edge set of LG^/A„(a)o is 

(2.6) : jGZ/nZ, /i=l,...,A;}, 

where the length of any edge (J~^(j), J~^(j + a/i)) is Ih- Hence J yields an isomorphism 
between the digraphs LG^ / Kn{a)Q and C^(£, a), preserving the quasimetric. □ 

2.2. Undirected lattice graphs. The discussion of the previous section applies with very 
small changes to the undirected lattice graph LGk = (Z'^ , E) , where the edge set E is the 
same as before but the edges are considered without orientation. 

The metric on LG^ is defined as for LG^, and now the distance between vertices m, n £ 
LGk is given by 

(2.7) (i(m, n) = (n — m)+ • 

where we denote Zj^ := {\zi\, . . . , l^^j) for any z = {zi, . . . , z^) G M'^. Furthermore if A is a 
sublattice of Z'^ then the distance between vertices m + A and n + A in LG^/ K is given by 

(2.8) d{m + A, n + A) = min{z+ ■ i : zGn-m + A}. 

Now take (a, n) = (ai, . . . , a^, n) G N*^ with < ai < . . . < Cfc < f , and recall the definitions 



(2.4) and (2.5). 



Lemma 2. The metric graphs LGfe/A„(a)o and Cn{i,a) are isomorphic. 
The proof is the same as for Lemma [T} 

2.3. Diameters. Let A be a sublattice of Z'^ of full rank (viz., of finite index). In view of the 
definition of the distance on LGj^/A, we have for the diameter 

(2.9) diam(LGfc/A) = max min|z+-£ : zGm + Aj. 

meZ^/A 

We define a corresponding diameter for the continuous torus M'^/A: 

(2.10) diam£(M''7A) = sup min{z+ ■ £ : z e y + A} . 

j/eiRfc/A 

This is the maximal distance between any two points on M'^/A, when distance is measured 
in the "f-weighted ^-"^-metric" , i.e. we define the distance between any two points a; + A and 
y + A on M^/A as the minimum of z^ ■ £ taken over all z £ y — x + A. 
Similarly for the directed graph LG^/A we have 

(2.11) diam(LG+/A) = max min ((m + A) n Z^g) • i- 

meZ'=/A 

We define a corresponding directed diameter for the continuous torus M.^ /A: 

(2.12) diamKM'^/A) = sup min ((y + A) n M|o) • ^• 



This is the maximal distance between any two points on M'^/A, when distance is measured in 
the ^-weighted £^-metric, and we only allow paths with non-negative components. 
Recah that we write e = (1, . . . , 1) G R^. 

Lemma 3. Let (a, n) = (oi, . . . ,ak,n) £ N'^ with < ai < ■ ■ ■ < < n. Then 

(2.13) diam(LG+/A„(a)o) = diam+ (M'=/A„(a)o) - e • 
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If furthermore fflfc < ^ then 

(2.14) diamf(MVA„,(a)o) - ^ < diam(LGfc/A„(a)o) < diam^(MVAn(a)o). 

Proof Set A = A„(a)o. Let y e R'^ he arbitrary. Set m := ([j/ij, . . . , [yk\) G Z'^, so that 
y = m + z for some vector z € [0, 1)*^. Using A C we have 

(2.15) (y + A)nM|o = -2 + (('^ + A)nZ|o), 
and thus 

(2.16) mm{{y + A) n M^q) • ^ = z • ^ + min((m + A) n Z|o) • £. 

Taking the supremum over all y G M.^, or equivalently the supremum over all {m,z) G 
iJ' X [O,!)'', we obtain 

(2.17) diam+(]RVA) = sup z • £ + diam(LG+/A) = e • £ + diam(LG+/A), 

z6[0,l)'= 

and we have proved ([2.13^. 



We next turn to (2.14). The right inequality in (2.14) is obvious from (2.9) and (2.10). 
To prove the left inequality, let y = (yi, . . . ,yk) be an arbitrary point in M''' . Then there is 
an integer vector m = (mi, . . . , m^) satisfying \mj — yj] < ^ for j = 1, . . . ,k. Now for any 
z £ m + A there is a point z' G y + A satisfying | z'j — Zj\ < ^ for all j . Hence 

(2.18) 

e ■ i e ■ £ 

min{z+ ■ £ : z £ y + A] < min{2+ ■ i : z £ m + A] -\ < diam(LGjt/A) H . 



Since this holds for all y G M'^ we obtain the left inequality in (2.14). □ 
Now set 

(2.19) L>„(£) = diag(n-i/''4,--.,n-i/^'4) G GL(A;,M), with U = nh---£k. 
We have det-D„(£) = n~^, and hence 

(2.20) L„^a,i ■■= Ania)oDn{£) 



is a lattice in M of covolume one, viz. Ln^a,e. £ ^k- It is also clear from the definition (2.10) 
that this transformation translates diam^(M'^/A„(a)o) into an unweighted (or "e- weighted" ) 
f^-diameter, viz. 

(2.21) diam^ (MVA„(a)o) = diame (MV^n.a/) • 
Similarly 

(2.22) diam+(M''7A„(a)o) = diam+ {R''/Ln,a,e) ■ 
Combining Lemma [T| Lemma [2] and Lemma [3j we have now proved: 

Proposition 1. Let (a, n) = (ai, . . . , ak,n) G N*^ with < oi < • • • < < n. Then 

(2.23) diamC+(£, a) = m diamt {^^ / Ln,a/) -e £. 
If furthermore < ^ then 

(2.24) diame (MVin,a,£) - ^ < diamC„(£,a) < m diam^ (MV-^^n.a/) • 
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2.4. Diameters and covering radii. We next note that, for an arbitrary /c-dimensional 
lattice A C M^, the ^"'^-diameters diame(M^/A) and dianig (M'^/A) can be interpreted as the 
covering radius with respect to A of the simplex A and the cross-polytope respectively. 
(Recah ([TS]) and ([TTs]).) 

Lemma 4. For any lattice A C M'^ of full rank we have 

(2.25) diame(M'=/A) = p{^, A) 
and 

(2.26) diam+(RVA) = p(A, A). 
Proof. Note that, for any y € M^, 

(2.27) min{z+ • e : z £ y + A} = sup{R > : n (y + A) = 0}. 

Hence by (2.10), diame(M'^/A) equals the supremum of all -R > such that there exists a 
translate of A which is disjoint from One sees that this holds if and only if + A 7^ M'^. 
Hence 

(2.28) diame(M'=/A) = sup{R > : R^ + A^R^}= p{^, A). 



The proof of (2.26) is the same, using the fact that 

(2.29) mm{{y + A) n M|o) • e = sup{i? > : i?A n (y + A) = 0} 

for ah y gR''. □ 

2.5. Equidistribution. The key to the proof of Theorems [T] and [2] is the following equidis- 
tribution theorem, which is a consequence of Theorem 7 in [28] . 

Theorem 3. Let I = (^1, ... ,4) S R^q, and let V d {x G R'^ : ^ < xi, . . . ,Xd-i < Xd} be 
a bounded subset with boundary of Lebesgue measure zero. Then for any bounded continuous 
function f : — )• R, 

In order to prove Theorem [3] we first prove Theorem |4] below, which is a corollary of [28^ 
Thm. 7]. Set G = SL{d,R), Gq = SL(fc,M) and T = SL{d,Z), Tq = SL{k,Z). For any M G Go, 
Tj^M is a /c-dimensional lattice of covolume one in R^, and this gives an identification of 
with the homogeneous space ro\Go. Then is identified with the unique Go-right invariant 
probability measure on ro\Go; we also use the same notation /xo for the corresponding Haar 
measure on Gq. Let H be the following subgroup of G: 

(2.31) H = \m = ^ ) : A G Go, c G 1 



We normalize the Haar measure fiff of H so that it becomes a probability measure on T\TH; 
explicitly: 

(2.32, ..„(M)^.„(.,.c. M^(^ f 

where dc denotes the standard Lebesgue measure on M^. We set 

(2.33) D'M--= (^0^^^ ^) =diag{U~y%,...,U~'/'ek,n)eG iU = nh- ■■£,). 
Theorem 4. 

(i) For every {a,n) G N'^ we have u{n^^a)D'^{£) G TH. 
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(ii) For any I G M^,q, any hounded subset V C {x £ : < xi, . . . , Xd~i < Xd} with 
boundary of Lebesgue measure zero, and any bounded continuous function /o : T\TH — )• 
M, we have 

(2.34) ^lim fo{u{n-'a)D'M) = h{M) d^M). 



Proof. To prove (i), note that for any (a,n) G N"^ there exists 7 G T such that {a,n)-f = 
and for this 7 we have 



(2.35) V(?^"^a) 



A H) 

c n^^ 



for some c G M'' and A G GL(/!;,M) with det A = n. It follows that ^^u{n^^a)D'j^{£) G H, and 
hence it(n-ia)D;(£) G TH, 

Next to prove (ii), note that since T\TH is an embedded submanifold of T\G, it suffices 
to prove that ( 2.34[ ) holds when /q is an arbitrary bounded continuous real- valued function 
on T\G. But this follows by applying Theorem 7 in [28^ with the test function f{x,M) = 
fo{MD'^^{i)). □ 

Proof of Theorem We have 

(2.36) Ln,a,e = {{^^ X nZ)n(a))pZ)„(£) = {z''+'u{n-'a)) ^Dn{£) = {lI'+^u{n-^a)D'^{t))^. 
Hence the left hand side of ( 2.30| ) can be expressed as 

(2.37) ^lim^ f[{^'^Mn-'a)D[Xt-)),). 

{a,n)eN''nTV 

Let us now apply Theorem |4] with the test function /o given by fo{M) := /((Z'^"''^M)o) for 
all M G TH. To see that this is well-defined, note that ifM = 7(^f) with 7 G F and 
( ^ *0 ) G then 

(2.38) (Z'=+^M)o = (^^ = Z''A G Xk 

so that fo{M) = f{'L^A). The function /o is obviously bounded and F-left invariant; further- 
more the formula fo{M) = fifi' A) just proved shows that /o is continuous on Tl, and hence 
on TR. Now Theorem |4] gives that the limit in (|2.37|) equals 

= /(Z'A)d„(A), 

and we are done. □ 

Theorems [T] and [2] now follow from Theorem [3] combined with ( |1.3 ), Proposition [T] and 
Lemma |4j Indeed, let £ G M^q and P C 5"*" be given as in Theorem l| Then Theorem [3] 
together with ( |1.3[ ) imply that if we view L'^i ^ ^ as a- (Xfc-valued) random variable defined 
by taking (a, n) uniformly at random in N'^'"*""'^ n TT>, then as T — )• 00, Ln^a,£ converges in 
distribution to a random variable L G Xf. taken according to hq. We next note that the 
functions L 1— )• p(^, L) and L 1— )• p(A, L) are continuous on (this is immediate from [16^ 
Prop. 4.4]; for the case of A it was also proved in |28[ Lem. 4, Thm. 9]). Hence by the 
continuous mapping theorem, 

p{^,Ln,a,£) piV,L) and p{A, Ln,a,i) pi^, L) asT-^oo. 
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Thus by Lemma 4, and using the obvious fact that {n£i • • • £k) ^ -— 0, we have both 
diam+(MVW) - (^^^.".^^)v^ ^ P(A,i), 

diame(MV^n,a,^) " ^, / " A p{^, L), and diame(MV^n,a/) A L) 

as T — )• oo. Hence Theorem [T] follows from Proposition [l| and so does Theorem [2] if we also 
assume V C ^. □ □ 

3. On THE DISTRIBUTION OF L) FOR RANDOM L £ Xk 

In this section we give proofs of those results about the distribution of p{^, L) for random 
L € Xk which we have used or mentioned in previous sections. 

3.1. Proof of the continuity of Pk{R). The proof that Pk{R) is a continuous function of 
R follows the same basic strategy as the proof of the continuity of Pk{R) = ^d{R) in [28^ 
Lem. 7], but the details are a bit more complicated. We start by giving a necessary criterion 
for p(^,L) = R, in Lemma [5] below. We write {±1}'^ for the set of all vectors in M'^ of the 
form ibei it 62 ± ... it e^. For each e = (ei, . . . , e^) E {±1}'^ we let be the (closed) face of 
*P given by 

(3.1) ^^ = {xe^i : e-x = l} 

k 

= |a; = (xi, . . . , Xk) G : ejXj = 1 and ejXj > for all j = 1, . . . , A;|. 

i=i 

It is clear from the last relation that is a (A; — l)-dimensional simplex. The faces Cpe 
together cover the boundary of 

k 

(3.2) d^} = {xe M'- : J2 l^il = = U 

i=i ee{±i}'= 

Lemma 5. // /o(*P, L) = R for some L G Xk and R > 0, then there is a vector G M'^ and a 
nonempty subset E C {±1}'^ such that 

(i) L n (C + i???°) = 0; 

(ii) L n (C + i?*Pe) / /or eac/i e G 

(iii) t/iere does not exist any q G M*' satisfying e • o; < for all e ^ E. 

Proof. Let L G and ii > be given with /?(*P, L) = R. Then, similarly to what we noted 
in the proof of Lemma [4| i? is the supremum of all R' > such that there exists a translate 
of R'^i which is disjoint from L. Hence by a simple compactness argument there is some 
translate of R'^° which is disjoint from L, i.e. we have L n {(^ + R'^°) = for some ^ G M'^. 
Let us fix such a vector C, and let E be the set of all e G {il}'^ for which Ln{C + R!^e) 7^ 0- 
Then conditions (i) and (ii) hold by construction. Assume that (iii) does not hold, and let a 
be a vector in M'^ satisfying e • a < for all e G -E. Now because of L fl (C + R^°) = and 



(3.2) , for every point a; G (C + -R*P) H L there exists some e £ E such that a; G C + R^^e- In 
particular we then have e ■ {x — = R, and hence, for all t > 0, 

(3.3) e ■ {x - (C + ta)) = R-te- a > R, 

so that X ^ {C + tcx) + R^. It follows that (C + ta) + R^ is disjoint from (C + R^) n L, for 
every t > 0. Hence for t > sufficiently small, {C + ta) + R^ is in fact disjoint from all L, 
so that L n ((C + ta) + ii'q?) = even holds for some R' > R. This gives /j(*P,L) > -R, a 
contradiction. Hence also condition (iii) must hold. □ 

Lemma 6. A finite nonempty subset E = {ei, . . . , e^} C M*^ \ {0} satisfies the condition (iii) 
in Lemma^ if and only if Yll=i ^i^i = holds for some choice of ci, . . . ,Cr > 0, not all 0. 
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Proof. Let C be the conic hull of —E. Then condition (iii) in Lemma [5] says that the dual 
cone C* has empty interior, or in other words that C* is contained in a proper linear subspace 
of M^. This holds if and only if C** = C contains a line through the origin, i.e. if and only if 
Yll=i ^i^i — holds for some choice of ci, . . . , Cr > 0, not all 0. □ 

The following lemma shows that Pk{R) is continuous. 
Lemma 7. For every R > 0, 

(3.4) fio{{L e Xk : p{V, L) = R}) = 0. 

Proof. By the definition of /xq, it is equivalent to prove that the set S of all A ^ Gq satisfying 
p{^,1j^A) = R satisfies po{S) = 0. Let £ be the family of all subsets E C {il}'^ satisfying 
condition (iii) in Lemma [5j Then, by that lemma, 5" is a subset of 

(3.5) [j{AeGo : there exists C G M'' such that Z^A n (C + R^e) / 0, Ve G ^}. 

But £ is finite; hence it suffices to prove that each individual set in the above union has 
measure zero. Thus fix some E £ £; say E = {ei, . . . , e^}. The corresponding set in the above 
union can be expressed as 

(3.6) y {AeGo : there exists C G K'' such that UiA G C + R'^e, (i = 1, . . . , r)}. 

ni,...,n,.GZ''" 

This is a countable union, and hence it suffices to prove that each individual set in the union 
has measure zero. Thus we fix some rii, . . . ,nr G Z^. Since E satisfies condition (iii) in 
Lemma [sj there exist, by Lemma [6| some ci, . . . ,Cr > 0, not all zero, so that Yll=i ^i^i — ^■ 
Now riiA G C + R^ei implies {riiA — ■ ei = R, and multiplying this relation with Cj and 
adding over all i we obtain X^^^^ CiUiA ■ = RYl'^i=i Cj- Hence the set corresponding to our 
fixed ni, . . . , n.r in the above union is a subset of: 

r r r 

(3.7) {^GGo : ^Cini^-ei = i?^Ci} = {^GGo : tr(MA) = i?^ Ci}, 

i=l i=l i=l 

where M = {m£j) is the k x fc-matrix given by m£j = Yll=iCi{ni ■ ej)(ej • eg). We have 
X]i=i Cj > 0, since ci, . . . , Cr > and at least one Cj is positive. Hence if M = then the set 



(3.7) is empty. If M 7^ then the set (3.7) is a submanifold of Go of codimension one (cf. the 
proof of [28„ Lem. 7]). Hence the set (3.7) has measure zero also in this case and the proof is 
complete. □ 

3.2. Proof of pk > ^(fc!)^/^ for A; > 3. We noted in {[TTsl) that pk > vol(«p)~i/^ = ^(A;!)^/'^ 



and in the present section we will prove that strict inequality holds in this relation when A: > 3. 



Since the infimum in (1.17) is known to be attained (cf., e.g., [211 Thm. 21.3]), it suffices to 
prove that there does not exist any lattice covering of M'^ by translates of ^ which has density 
exactly one, viz. with the *P-translates having pairwise disjoint interiors. In fact we will prove 
the stronger fact that there does not exist any tessellation (lattice or non- lattice) of by 
translates of ^: 

Proposition 2. For k > 3 there does not exist any subset P C M*^ such that P + = M*^ and 
(r + «p°) n (s + ^3°) = for allr ^ s £ P. Hence in particular, pk > ^(A;!)^/'' for k>3. 

The proof of this fact is quite easy but we have not been able to find an appropriate 
reference for it. The question of finding the optimal lattice covering of M'^ by translates of 
*P was studied by Dougherty and Faber in [16l Sec. 7], and they conjecture that the optimal 
density is |, which would mean that p3 = |iJ/2 = 0.9449 .... We remark that for the more 
classical question of lattice sphere coverings, the optimal coverings are known in dimensions 
up to 5; cf. [ITj, I26J, t40j. 
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Proof of Proposition^ Assume P + <p = M*^ and (r + ^3°) n (s + = for all r 7^ s G P. 
Without loss of generality we assume € P. Now for any point x = {xi, . . . , Xk) € M^q with 
xi + . . . + Xfc = 1 (i.e. X G we may argue as follows. For any e > we have x + ee ^ 
and thus x + ee £ r + ^ for some r = (ri, . . . , r^) G P \ {0}. Letting e — )• it follows that 
there exists a point r = (ri, . . . , r^) G P \ {0} such that x £ r + i.e. Ylj=i ~ — 
and also a; + ee S r + *p for all e's in some sequence of positive numbers tending to 0. Now 

k k k 

(3.8) I^^-I ^ E 1^^- - ^J-l + E^^- = 1 + 1 = 2' 

j=l j=l J,=l 

and if rj < Xj would hold for some j then we would have strict inequality in the above 
computation, and this would lead to the contradiction G n (r + To sum up, we 



have proved that for any given x G MJI^q with xi + . . . + x^ = 1, there exists some r G P \ {0} 



satisfying rj > Xj for j = 1, . . . , k, and X]j=i = 2. 

Let us first apply the above fact with x = (1 — (k — l)e, e, . . . ,e) with e > tending to 
zero. It follows that there exists some r £ P with ri > 1, r2, . . . , > and Yl'j=i ''"j — 2- 
Next we apply the above fact with x = (ri — 1 + e, r2 + e, rs — 2e, r4, . . . , r^) (here we use 
k > 3!). This leads to the conclusion that there exists some s G P with si > ri — 1, 
S2 > r2, Sj > rj for J = 3, . . . , k, and J2j=i — 2- In particular r ^ s since S2 > ^2. Now 
si = 2 — X]j=2 '^i < 2 ~ Si=2 '"i ~ '^i' ^^"^ hence 

k k 
(3.9) ^ \sj - rj\ = (n - si) + ^{sj - rj) = 2ri - 2si < 2ri - 2(ri - 1) = 2, 

i=i j=2 
which leads to the contradiction ^(r + s) G (r + *p°) n (s + q}°). □ 

3.3. The asymptotic formula for Pk{R). We now discuss the proof of the asymptotic 
formula stated in Remark 11.71 viz. 



(3.10) PkiR) = ^)+Ok{R~'~'~^-). 

It turns out that most of the proof in [38] of the asymptotic formula for Pk{R), carries 
over with very small changes to the present case: Mimicking [381 Sec. 2.1-3] we obtain 

(3.11) Pm = ^ f + o.(fl-"-«'- A), 

where S^^""^ is the unit sphere in M*^ centered at zero, dv is the {k — l)-dimensional volume 
measure on S^~^, and £{v) is the width of *p in the direction v, viz., for t; = (t;!, . . . , Vd) G S^^"^, 

(3.12) i{v) = 2ma^{\vi\,...,\vk\). 



Now to get (3.10) it only remains to prove the following. 
Lemma 8. For every k >2 we have 

(3.13) / e{v)-''dv = k. 
Proof. Let ^* be the polar body of i.e. 

(3.14) ^* = {x€R'' : x-y <l,yy e"^} = {rv : v£ S^-\ < r < {li{v))-^}. 
Then clearly 

(3.15) / £{v)-^dv = 2-^kvol{^*). 

However one verifies easily that equals the A:-dimensional cube [—1,1]^^; hence vol(*p*) = 2^^ 
and the lemma follows. □ 
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One may note that for k = 2, (3.10) says P2{R) = -^R ^ + 0{R which is consistent 



with the exphcit formula stated in Remark] 1.8 



4. The explicit formulas for P2iR) and p2iR) 



We now prove the explicit formula for the density P2{R) which we stated in Remark 1.8 
Proposition 3. For k = 2 the density Pk{R) = ~jiiRk{R) is given by 

_/0 (0<fl<-^) 

4.1. Auxiliary lemmas. To prepare for the proof of Proposition [3] we first prove a series of 
lemmas. As a first step, note that since ^ for A; = 2 is a square with side \/2, the formula 
( 1.15[ ) may be rewritten as (using the S0(2)-invariance of po) 



(4.2) P2{R) = no{{L e X2 : p{K,L) > r}), where r := V2R, 
and where K is the unit square 

(4.3) K := [0, if = {x = ixi,X2) : xi, G [0, 1]}. 

We will make frequent use of the fact that, just as in proof of Lemma |4| p{K,L) is the 
supremum of all r > for which there exists a translate of rK that is disjoint from L. 

We now introduce a parametrization of X2 that is tailored to give a practicable expression 
for (4.2). To motivate our definition below, note that by Lemma [5] (transformed from *P to K), 



if L G X2 and p{K, L) = r then there is some ^ G such that L has no point in the interior 
of + rK, but L has a point on each of two opposite sides of C + rK. By perturbing <^ in a 
direction parallel to these sides we may also, at least for generic L, assume that L intersects 
one more side of ^ + rK. If we assume that the three L-points on the sides of ^ + rK are 
C + r(a, 0), C + '^(O, P) and C + ^(1) 7) with a, /3, 7 G (0, 1) then it follows that L contains the 
vectors r{a, —f3) and r(l,7 — /3), and in fact these two vectors necessarily span L, since L is 
disjoint from the interior of + rK. Using also the fact that L has co-area one, it follows that 

(4.4) L = L(,,^,^) := <5-i (Z(a, -/3) + Z(l, 7 - /?)) 

where 

a -/3 



(4.5) d = (5(a,/3,7) : = 



1 7-/3 



(1 - a)(3 + a7 > 0. 



Lemma 9. The map (a,/3,7) 1— )• -Zj(Q,,^,-y) is a local diffeomorphism from (0, 1)^ to X2, under 
which the measure po corresponds to 

(4.6) ij5(a,P,-f)-^ dad/3 d-f. 



Proof. Set 



1 



(4.7) .^A.„,,,:.r.(t ^-/^).G„, 

SO that L(-^ 13 ■y) — A. A computation shows that the Iwasawa decomposition of A is given 

by 



(4.8) 
where 



1 x\ ( ^ \ ( ^ ~ ^^'^ ' 
1 M II Jy) Uin0 cos(/> 



, . a-B'y + B'^ (l-a)/3 + Q7 , vr x 

(4.9) - = YTW^--- ^=-iTW^'- 2 +=>-<-("-* 
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One furthermore computes 

(4.10, tff-<^-<^-)^)- 



It is clear from ( |4.9[ ) that x,y,(f) are smooth functions of {a,/3,j) E (0,1) , and since also 
the Jacobian determinant (4.10) is non- vanishing for all these (a,/3,7) it follows that the 
map (a,/3,7) i— )• {x,y,(j)) is a local diffeomorphism from (0, 1)^ to M x M>o x (0, vr). However 
the Iwasawa decomposition is known to be a diffeomorphism from (x, y,4>) G M x M>o ^ 
(M/27rZ) onto Go, under which the measure corresponds to dxdydcj). Hence the map 

(a, /3,7) 74(Q,^^ ,y-) is a local diffeomorphism from (0, 1)'^ to Go, under which |Uo corresponds 
to (4.6). To complete the proof of the lemma we need only recall that the quotient map 
Go — >■ X2 = ro\Go is a local diffeomorphism and fj,o on X2 is just the measure corresponding 
to fiQ on Go. □ 

Set 

(4.11) := Z(a, -/3) + Z(l, 7 - /3) C 

so that = 5~^L'^^ ^ By construction the translated lattice (0, /3) + L'^^ ^ contains 

three points on the boundary of the unit square K, namely (a, 0), (0, /3) and (1,7). We next 
determine those (a, /?, 7) for which (0, /3) + L'^^ ^ contains no other point in K. 

Lemma 10. Given (a,/3, 7) € (0,1)^, the relation 

(4.12) ((0, /3) + ,^)) n = { (a, 0), (0, /3), (1,7)} 
holds if and only i/ /? + 7 > 1 . 

Proof. We have 

(4.13) (0,/3) + L'(„^^ ,^) = {(0,/3) + ni (a, -/3) + 712(1,7-/3) : n = (m, na) G Z^}. 

In this representation the three points (a,0), (0,/3), (1,7) correspond to n = (1,0), n = 
(0,0) and n = (0,1), respectively. Taking n = (—1,1) in ( |4.13 ) we see that (4.12) implies 
(1 — a,/3 + 7) ^ -ftT, viz. /3 + 7 > 1. Conversely, assume /3 + 7 > 1. One then immediately 
checks that, in the above representation, those n with \ni\ < 1 which give points in K are 
n = (1, 0), (0, 0), (0, 1), and no others. To conclude the proof of the lemma it now suffices to 
show that all n G with |ni| > 2 also give points outside K. Assume the opposite, i.e. that 

(4.14) p := (0, (3) + ni(Q, -/3) + n2(l, 7 - /3) G 

for some n G Z^ with |ni| > 2. It is a simple geometric fact that for any such n, there exists an 
integer m such that the point (sgn(ni),m) belongs to the closed triangle with vertices (0,0), 
(0, 1) and n. Applying the affine map {x, y) 1— )• (0, /?) + x(a, — /3) + y(l, 7 — /3) we conclude 
that the point 

(4.15) q := (0, /3) + sgn(ni)(a, -/3) + m(l, 7 - /3) 

lies in the closed triangle with vertices (0,/3), (1,7), P- Hence, since q is not equal to one of 
the triangle vertices, and since < /3,7 < 1 and p £ K, we conclude that q G K°. This is a 
contradiction since we saw above that no point in (4.13) with |ni| < 1 lies in K°. □ 

Set 

(4.16) 0:= {(a,/3,7) G (0,1)3 . ^ + 



After a translation and a scaling, Lemma 10 says that for any (a, /3, 7) G ri, the lattice L((^^^ 
meets 6~^{0, — /3) + 6~^K in exactly three points, all lying on the boundary of this square. 
Hence for such (a,/3,7) we have p{K^L(^^^p^^^) > (^~2. The next lemma shows that we always 
have equality in this relation. 
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Lemma 11. For any (a,/3,7) £ 0, we have L'^a i3 -y) ^ iC + rK°) / for all r > I, C, £ R^, 
and hence p{K, Lf^^^f^^^^) = 6~2 . 

Proof. Assume the contrary; then Lj^ ^ H + rK°) = for some r > 1 and ^ G M^. Note 
that there exists t > such that L'^^ ^ H (C — (0, t) + rK°) ^ (this follows since Lj^ ^ 
contains a vector with positive ei-component < r, e.g. the vector (a,0)). Taking to ^ to be 
the infimum of all t > with that property, and then replacing ^ with ^ — (0,to), we obtain 
a situation where the side {C + {x,0) : < x < r} contains a lattice point i £ L'^^ ^ while 
stm L[^^^^^^ n (C + rK°) = 0. But now also £ + (-a, /3) G L[^^^^^^ and £ + (1 - a, 7) g' Lj^^^ , 
and at least one of these two points must lie in + rK°, since £ £ {C + ix,0) : < x < r} 
and r > 1. This is a contradiction. □ 

Lemma 12. T/ie map (a,/3,7) 1— )• L(a^^ .y) is a dijfeomorphism from 0, onto an open subset 

X'2 0fX2. 



Proof. In view of LemmalOlit suffices to prove that the map is injective. Thus assume ,),) 
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-^(o',/3',7') for some (a, /3, 7), (a', /?', 7') G ^2. Then (5(a,/3,7) = 6{a' , /3' ,^') by Lemma 
hence L'^^^^^ = L'^^, ^, Call this lattice L'. Using now Lemma 10 and (a, — pj G L' 
it follows that (a,0) + L' is disjoint from K°. In particular (a,0) + [—a',(3') ^ and 
(a, 0) + (1 — a', 7') ^ and these two relations together imply a' = a. Now also f3' = f3 and 
7' = 7 follow easily. □ 

Let W = (10); this element acts on by switching coordinates, and it acts on lattices 
L C by L h-^ LW := {xW : x € L}. The latter action gives a diffeomorphism of X2 onto 
itself, preserving fiQ. Set 

(4.17) X'^ := X'^W, 
where X2 is the open subset of X2 defined in Lemma 
Lemma 13. X'^ n X'{ = 0. 

Proof. Assume the contrary; then L(^a,i3,-y) = some (a, /?, 7), (a', /3', 7') G O. 

Now p{K, L(^^i^iji^y^W) = p{K, L(^^i^f^i^y^), since W maps -fC onto itself; hence by Lemma 11 
we have 5{a, /3, 7) = 6{a' , /?', 7'), and thus also L'^^ fB ^) ~ ^'{a' /?' 7')^' ^^^^ ^^^^ lattice L'. By 
Lemma fTOl we have 

(4.18) ((0, P) + L')r^K = {(a, 0), (0, /3), (1,7)}. 

Using here (7', 1 — q') G V we get (0,/3) + (7', 1 — a') ^ K, viz. j3 > a'. On the other hand 
using (/3', — a') G L' we get that (0, /3) + (/3', — q') is either outside K or else equals (a,0); 
hence we must have /3 < a'. This is a contradiction. □ 

Lemma 14. /io(X2 U X'^) = 1. 
Proof. We have 



(4.19) po{X'2UX'^) = 2po{X'2) = ^ 6-'dadpdj- 

vr^ 7n vr^ 7n ((1 - a)p + a7)' 

by Lemma [9j Writing this as an iterated integral and evaluating the innermost integral over 
a G (0, 1), we get 

(4.20) 

6 d/3d-/ _ 6 /"^ log(l - 7) ^_ _ 6 p X 6 

'1-7 



211 R - 21 d^=^ -— dx = ^r{2)C{2) = 1. 

vr^ Jq Ji f3j vr^ 7 vr^ Jo - 1 vr^ 



(We substituted 7 = 1-6"^ and then used [22l Thm. 14].) □ 
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Remark 4.1. Another way to prove Lemma 14 is to make the discussion preceding (4.4) more 
precise, so as to show that for a generic lattice L G X2, there exists some ^ G such that 
LCi {C + rK°) = and either L contains the three points C + r{a, 0), C + ^^(0, /3) and C + r{l, 7) 
for some {a,f3,j) G i7 (thus L G X2), or L contains the three points C, + r(0, a), C, + r{l3,0) 
and C + '^(7) 1) for some (a, /3,7) G (in which case L G X2). However the above proof by 
direct computation also serves as a nice consistency check of our set-up. 



4.2. Proof of Proposition 3. Using (|4.2|), Lemmas |9[ ll| 12| 13, 14 and the fact that 

da 



L I— LW preserves both fiQ and p{K,L), we get: 



(4.21) 



P2{R) 



6 



where r = y/2R, 
(4.22) 
and 



6 



1 rl 



7r2 Jq Ji_^Jj^ ,^^^ (5(a,/3,7)2 
nr := {(a,/3,7) GO : 5(a, /3, 7) < r'^}, 
//3,7,r- = {a e (0,1) : (5(q,/3,7) < r-2|_ 



d(3 d'-f, 



(4.23) 

Recalling (5(a, /3,7) = (1 — a)/3 + 07 we find that //3,7,r = (0, 



if ^ < r-2 < 7, 



,7,r 



(- 



-/3 



7-/3 



, 1) if 7 < r ^ < /3, while 



/3,7,r 



(0,1) if /3,7 < r-2 and 1/3, 



,7,r 



if /3,7 > r ^. Now 



it is easy to compute the derivative of the innermost integral in (4.21 ) with respect to r, using 



the fact that 5(a,/3,7) 
(4.24) 



dr 



for a 
da 



. If /3 < r < 7 then we get 
d r-2 _ , 2r 



/3,7, 



, (J(a,/3,7)2 Vdr 7-/3 



7-/3 



Similarly when 7 < r ^ < /3 we get 



(4.25) 



d 
dr 



da 



2r 



,^^^,.<5(a,/3,7)2 /3-7' 

while if /3, 7 < or /3, 7 > then the derivative vanishes. Hence we obtain, using also the 
symmetry between /3 and 7: 

12^/2 



(4.26) 



2r 

7^ 



(i7. 



where Jr is the set of all pairs (/3,7) G (0, 1)^ satisfying both /3 + 7 > 1 and /3 < r^^ < 7. If 
r < 1 then = 0, so that P2(-R) = 0. On the other hand if r > \/2 then we get 



(4.27) 



P2iR) 



24^/2 







1 d7 d/3 24^/2 



7-/3 



log 







1-/3 
1 - 2/3 



d/3 



I2V2 



-^r((l - 2r"2) log(l - 2r"2) - 2(1 - r'^) log(l 
Finally if 1 < r < \/2 then we get 



(4.28) p2iR) 



24\/2 



24^/2 



l-r" 



1-r- 



log 



d-fdp 
1-/3 



+ 



l-r-2 

(i/3 + 



1 d7d/3 
-27-/3 

log 



1-/3 



„-2 



I2V2 



■ 1 — 2/3 / j]^_^-2 
f (1 - 2r~2) log(2r-2 - 1) - 2(1 - r''^) log(l 
Hence, recalling r = \piR^ we obtain the formula stated in (4.1). 



/3 



d/3 



□ 
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4.3. The explicit formula for p2{R). We next turn to the exphcit formula for P2{R) which 
we stated in (1.12). This formula is due to Ustinov [39], who proved it by an argument 
involving Kloosterman sums and continued fractions. We think it may be of interest to see 
an alternative derivation of (1.12) based on the definition of P2{R) in terms of Haar measure 
on the space of lattices, cf. (1.7), and so we give an outline of this argument here. 
The overall structure of the argument is similar to the previous case of p2{R)- 

For any (a,/3,7) G (-5, \f we set A(a,/3,7) := ^~^^(a,/3,-y)! where 
(4.29) A'(„^^^^) := Z(-i + ^, -a - 7) + Z(/3 + 7, - 

and 



7) 



(4.30) 



K = k(q;, /3, 7) := 



-2+7 
/3 + 7 



-a 
1 
'2 



i + + a7 + /37 > 0. 



The motivation of the above definition is that the translated lattice ~ 7) 5 + 7) + ^ 

constains the three points (0, ^ — a), (2 + /?) 0) and — 7, 2 + 7) on the boundary of A. 

By a similar computation as in Lemma [9] one proves that the map (a,/3,7) 1— )• A(q,^^ ,^) is a 
local diffeomorphism from 

(4.31) 



1 l^3 

21 2' 



to X2, under which the measure fiQ corresponds to 
;K{a, 13, 7)"^ da d(3 d^. 



Next one proves analogues of Lemma [To] and Lemma [TT] It is useful to assume that at 
least two of a,f3,^ are positive. Note that for generic (a,/3,7) G (~2) 5)'^ we can always get 
to this situation after possibly applying the map W = (5 J) (cf. Section 4.1); this is because 



AW = A. 



and A(„_^^^)H^ 



Ar 



-7)- 



It now turns out that if 



(a,/3,7) € (— |, |)^ and at least two of a,/3,7 are positive, then the necessary and sufhcient 
condition for (5 — 7, 5 + 7) + p 7) contain no other points in A than (0, | — q), (^ + /3, 0) 
and (i - 7, 2+7), is: 

(4.32) a + /3>0, a + 7>0, /3 + 7>0. 



(Note that, in the other direction, (4.32) implies that at least two of a,/?, 7 are positive.) 
Next, for any (a,/3,7) G (—5; 5)^ satisfying (4.32), the necessary and sufficient condition for 
^(a,/3,7) n (C + rA°) / to hold for ah r > 1, C G K^, is a + /3 + 7 < Set 

(4.33) n := {(a,/3,7) G (-5,5)^ : a + /3>0, q + 7>0, /3 + 7>0, a + /3 + 7<^|. 



holds for all (a, /3, 7) G 0,. 



It then follows from the last statements that p(A, A^^^^^^) = 1 

It now follows by similar arguments as in Lemma 12 and Lemma 13 that the map (a, /3, 7) 1— >• 
A(a,^,7) is injective when restricted Q, and hence gives a diffeomorphism from 0, onto an open 
subset X2 C Ar2, and furthermore that X2 is disjoint from X2 := X2W. Finally, it turns out 
that the union of X2 and X2 has full measure in X2: 

(4.34) fio{X'2UX'^) = l. 

(This can be proved either by a direct computation, cf. below, or else by proving that a generic 
lattice in X2 indeed must belong to either X2 or ^2 •) 
Using (1.7) and the above facts, it follows that 

6 f 

/ K{a, I3,'j)~'^ da dp dj, 



(4.35) 



P2{R) 



where now 
(4.36) 



R 



{(a,/3,7) G n : K(a,/3,7) < R'^}. 



We next introduce s = a + /3 + 7 and t 



+ + 7^ as new variables of integration in ( 4.35 ) . 



Note that < s < i for all (a,/3,7) G O.; also t > is^ by Cauchy's inequality. Conversely, for 
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given s S (0, h) and t > is^, the set of corresponding points (a,/3, 7) G is the circle with 



center s, s) and radius yt— ^s"^ in the plane {a + /3 + 7 = s}, and we may parametrize 
these points as 

(4.37) (a,/3,7) = |(s,s,s) + \Jt- ^s^ (^{cosuj)bi + (sinw)b2), w € M/27rZ, 

where 61, &2 is an arbitrary fixed orthonormal basis in the orthogonal complement of (1,1,1) 
in M.^. Now (a,/3,7) € 0, holds if and only if < s < ^ and max(a,/3,7) < s, and the latter 
condition is equivalent to (a, /3, 7) lying inside a certain equilateral triangle with side 2\/2s and 
center |(s, s, s) in the plane {a + /3 + 7 = s}. This triangle has inradius and circumradius 



2y |s; hence if yt — < y |s (viz., t < s^) then aZZ correspond to points in $7, while if 

< \J't ~ < (viz., < t < 3s^) then certain subintervals of u; G M/27rZ have 

to be removed, and the Lebesgue measure of those oj G M/27rZ which correspond to points in 
Q is 



(4.38) 27r - 6arctan(^y|s" - . 

Hence, using also K(a,/3,7) = \ + \s'^ — \t and I ^Q{^stw) I ~ 2vl' obtain: 



27r — 6 arctan^y |s ^yjt — 



(4.39) +/ / .1 ^ _ i,,2 -dtds 



where 

(4.40) J^^/j = (I s^s^) n (s2 + 1 _ 2i?-2^ 00) and J,,r = (s^ 35^) n (s^ + i - 2i?-^ 00) . 

In particular for R < -v/3 we have Is^r = (|s^,s^) and Jg^/j = (s^,3s^) for all s G (0, ^) 
and in this case P2{R) = 1, corresponding to the fact that the union of X2 and X2 has full 
measure in X2, cf. \AM\ . Next if ^/3 < i? < 2 then stiU Js,r = (s^3s2) for all s G (0, |), 
but now Isfi = (gs2,s2) only for s G (0, ^yjAR-^ - 1], while = (s^ + i - 2i?-2,s2) foj. 
s G [^\/4i2~2 _ 1^ Hence by differentiation we obtain 

(4.41) p2{R) = -±P,{R)='^-^ R^.4R-Us = -(^-V4:^A. 

dR ^ ' vr 7x^V 4R-2-i vr / 

Finally if i? > 2 then /s^/j = for all s, and J^^r = for s G (0, - AR-"^], and J^^r = 
(§2 + 1 - 2i?-^ 3^2) for s G [^Vl -4i?-2, \). Hence by differentiation, 

(4.42) p2(R)='^r R* ■ 4i?-3 . 27r - 6arctanf #^1 - 4ii-2 • s"^)) ds, 

and this is easily evaluated to yield the expression given in (1.12). Hence (1.12) holds for all 
R>0. □ 



5. Further results 



We conclude by discussing a number of natural extensions and variations of Theorems [T] 
and [2] They require only minor modifications in the proofs. 
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5.1. Non-constant lengths. We now admit lengths I = (^i, . . . that depend on n and 
a. Such a requhement may arise for instance when C„(a) or C^(a) is embedded in a metric 
space (M^, say), and the lengths t are induced by the actual distance in that metric space. 
To make a precise statement: Let I : [0, 1]'^ — t- M>q be continuous, and assume £{x) > for 
(Lehesgue-) almost all x G [0, l]'^. Then, for any bounded set V C 5"*" with nonempty interior 
and boundary of Lebesgue measure zero, we have convergence in distribution 

(5.1) ^ ^p{A,L) asT^oo, 

{nii{n~^a) ■ ■ ■ £k{n-^a)) ' 

where the random variable in the left-hand side is defined by taking (a, n) uniformly at random 
in N^~^^ nTV , and the random variable in the right-hand side is defined by taking L at random 
in Xk according to fiQ. The analogous statement holds in the undirected case. 

The limit distribution of Frobenius numbers proved in |28] can be viewed as a special case 
of the above result, obtained by taking £{x) = x. Indeed, for this choice of i we have 

(5.2) diamC^(£(n~^a), a) = diamC^(a, a) = 1 + n~^F{ai, . . . , a^, n), 

where F{ai, . . . , a^, n) denotes the Frobenius number of the k + 1 numbers ai, . . . , a^, re; cf. 
[10\ Lem. 3] or |3l Sec. 2]. Because of this relation, and since the Frobenius number is invariant 



under permutation of the arguments, \28\ Thm. 1] follows from (5.1) 



5.2. The distribution of distances. Besides the diameter it is natural to consider the 
distribution of the distance between two randomly chosen vertices i and j. The ath moment 
(for a G Z>i) of this distribution is 



(5.3) = ^J^d(i,j) 



where re is the number of vertices. If A is a sublattice of Z'^ of finite index then in view of the 
definition of the distance on the directed quotient lattice graph LG^/A, cf. ( |2.3[ ), we get 

(5.4) M.[LG+/A]= J Yl (min((m + A)nZ|o)-^)". 



Similarly for the undirected quotient graph LGk/ A, we get via (2.8), 

(5.5) M^LGfc/A] = ^ (min{z+.£ : zGm + A})" 

Following the same strategy as for the diameter one can show that under the same assumptions 
as in Theorem [T| 

where 

(5.7) ^'L(y) :=min((2/ + L)nMy -e. 

Note that the scaling factor is the same as for the diameter, the maximum value of the 



distribution of distances; this is a non-trivial fact. In fact joint convergence holds in (5.6) for 



all a > 1, and from this it is possible to conclude that the distribution of normalized distances 



(n£i vertices z,j picked uniformly at random in C+(£, a), converges in distribution, 

as T — )• oo, to the distribution of ^'L(y) for y picked at random in M'^/L according to the 
standard volume measure dy. The convergence here is in the space of probability measures on 
M>o, cf., e.g., |25[ Ch. 10], and the setting of the limit relation is the same as in Theorem[T] The 
limiting random probability measure on M>o obtained in this result satisfies many interesting 
and beautiful properties; we postpone a detailed discussion of these matters to a future paper. 
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The analogue of (5.6) in the undirected case is 



K« M„(C„(Aa)) f 

(n<....4)°/^ ^ Xvx ' " 

where 

(5.9) ^L{y) ■■= min{2+ • e : z €y + L]. 

5.3. Shortest cycles. The shortest cycle length (scl) of a circulant graph and its connection 
to the geometry of lattices is discussed in [11] . The length of the shortest cycle in a directed 
quotient lattice graph LG^/A is 

(5.10) scl[LG+/A] = min(A n \ {0}) • £ 

In the undirected case, there are trivial cycles which correspond to cycles in the covering 
lattice graph LG^- The shortest of these have 4 edges, and thus the girth of any quotient 
graph LGk/ A is at most 4. We will ignore such cycles and only consider those which do not 
lift to a cycle in -LG^, or in other words cycles which have non-zero homology when viewed 
as closed curves on the real torus M'^/A. With this convention, the shortest length of all 
non-trivial cycles in a quotient lattice graph LG^/A is given by 

(5.11) scl[LGfc/A] = min {m+ •£ : m G A\{0}}. 

Using the same method as for the diameter one can show that, under the same assumptions 
as in Theorem [T} 

(5.12) '^'^^^^f^"?! ^ n M|o \ {0}) • e as T ^ oo. 

The complementary distribution function of the limit distribution in this relation is 

(5.13) Pfc,sci(i?)=/io({iG^fc : i?AnL\{0} = 0}), 

since for any L G Xk we have min(L n M^g \ {0}) • e > if and only if i?A n L \ {0} = 0. The 



analogue of (5.12) in the undirected case is 



(5.14) ^ ■ as T^oo, 
and here the complementary distribution function of the limit distribution is 

(5.15) Pfc,scl(^)=/^o({iG^fc : fiq3nL\{O} = 0}). 
Comparison with j30t Thm. 2.1] shows that for k = 2 the limit distribution in the directed 



case, ( |5.12[ ), (5.13), is related to the distribution of angles of two-dimensional lattice points 



(including multiplicities) via the formula 

(5.16) P2,sci(^) = ^o,o(0,a) 

with a = E? /2. Formula (2.16) in [SD] shows therefore that the density of P2,sc\{R) is related 
to the gap distribution function Po{s) for angles of lattice points, 

(5.17) P2,sci(i?) := -^P2Mi.R) = RPo{R^/2). 

An explicit formula for Po{s) can be derived from |4J (use Eq. (2.31) in |30] to relate Po{s) to 
Po{s); the latter is denoted Gd(s) in [1]); we find 



(5.18) Po{s) 



6 



1 (0<s<^) 
s-i(l + log2s) - 1 {\<s<2) 
s-i - 1 + Vl - 2s-i - 2s-^ log (^(l + \/l - 2s-i)) (s>2). 
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(5.19) P2,scl(i?) 



6 



Thus 

'R (0<i?<l) 

2R-^{l + 2\ogR) - R {1<R<2) 

2R-^ -R+ VR^ - 4 - 4i?-i log {l{l + Vl- 4ii-2)) {R>2). 

Similarly, [30, Thm. 3.1] shows that for k = 2 the limit distribution in the undirected case, 
( 5.14| ), (5.15), is related to the distribution of disks in random directions via the formula 



(5.20) /^2,.ci(i?) = i^o,o(0,a) 

with 2a = R?. To see this, note that (in view of the S0(2) invariance of hq) the square *P 
can be replaced by the square [-R/y/2,R/yj2f which in turn (due to the invariance under 
the symmetry x H> —x) can be replaced by the rectangle [0,i?/\/2] x [—R/\/2,R/y/2]. The 
function Fo,o(0, o") is in turn related to the free path length $o(C) of the two-dimensional 
periodic Lorentz gas via formula (4.3) in [30]. This implies for the density of P2,&c\{R)'- 

(5.21) P2,scl(^) := -J^P2,scliR) = RMRV"^)- 

The explicit formula for <I>o in [S] (denoted there by h; the formula can also be obtained from 
Eqs. (15) and (34)] or from [37, Prop. 3 ("r = 0")]) yields 



(5.22) P2,sci(ii) 



12 



R (0 < i? < 1) 

^(i+iog(2S)) a<R<V2) 

{R>V2). 
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